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Abstract 

State selective field ionization detection techniques in physics require a specific pro- 
gression through a complicated atomic state space to optimize state selectivity and 
overall efficiency. For large principle quantum number n, the theoretical models 
become computationally intractable and any results are often rendered irrelevant 
by small deviations from ideal experimental conditions, for example external elec- 
tromagnetic fields. Several different proposals for quantum information processing 
rely heavily upon the quality of these detectors. In this paper, we show a proof of 
principle that it is possible to optimize experimental field profiles in situ by running 
a genetic algorithm to control aspects of the experiment itself. A simple experiment 
produced novel results that are consistent with analyses of existing results. 



1 Introduction 



Many problems in experimental physics involve optimizing a set of param- 
eters in a large state space. Theoretical models exist for many problems of 
interest, but they are typically only valid with the application of a series of 
approximations, many of which do not necessarily hold true in a real labora- 
tory situation. Complications arise, for example, when a stray electromagnetic 
field is present that perturbs the system of interest. Such fields are often im- 
possible to remove completely from an experiment, and so it is frequently 
found that the measured parameters differ from the predictions. 



The field of quantum computing has expanded rapidly in recent years to en- 
compass many different fields with its promise of exponential or polynomial 
speedup in a number of important computations, such as factoring [T] and 
searching [2j. To build a useful quantum computer, it is necessary to gain 
exquisite control over individual quanta (e.g. atoms, electrons or photons) in 
a system. The advantage of quantum computing derives directly from the fact 
that extremely large state spaces can be manipulated in non-trivial ways to 
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perform certain calculations very much more efficiently than their classical 
counterparts. Evolutionary techniques have been successfully applied to some 
problems in quantum computing, in particular in the design of circuits to per- 
form quantum logic operations [3|. There exist several proposed schemes for 
implementing quantum computing, for example ion traps [I], nuclear magnetic 
resonance (NMR) [5] and cavity quantum electrodynamics [Bj. 



2 Quantum Computing in Cavity Quantum Electrodynamics 

The field of cavity quantum electrodynamics (cavity QED) has proved itself to 
be one of the leading candidates for quantum computing, since it offers many 
advantages, including the ability to couple light and matter in a controlled way. 
Current state of the art implementations of cavity QED experiments involve 
passing a sequence of carefully prepared atoms through a superconducting 
cavity that houses a microwave field containing small numbers of photons. 
For a number of technical reasons, so-called Rydberg states [7] are used, which 
means that a valence electron has been promoted via laser excitation to very 
high principal quantum number n. This leaves the electron loosely bound to 
the nucleus allowing the atom to be easily ionized in an electric field. The 
nature of the interaction is such that only two different Rydberg states partic- 
ipate, so the ability to distinguish between them after the cavity is sufficient 
to allow us to perform measurements for quantum computation. As in many 
other cutting edge atomic physics experiments, the atoms are typically alkali 
metals, in our case rubidium. 



3 State Selective Field Ionization 

Applying an electric field to an atom causes its energy levels to be Stark shifted 
to different energies. The operating principle of field ionization detectors is 
that it is possible to apply a varying electric field to 'walk' an energy level to 
the ionization limit, at which point the electron escapes and is collected in an 
electron multiplier, producing a detector 'click'. Different starting states reach 
ionization at different electric fields. Fig. [1] shows a simplified schematic of the 
effect of an electric field on the energy levels of the atom. Due to the interaction 
of the valence electron with the core of an alkali atom, features known as 
avoided crossings are exhibited in the Stark maps. The size of these avoided 
crossings defines the energy required to jump from one level to the next. If the 
electric field is ramped slowly with respect to this energy gap, the path will 
tend to follow an adiabatic path (dotted arrow in fig. [1]) and ionize exactly 
at the classical ionization limit. If the field is changed quickly with respect to 
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Fig. 1. Simplified schematic of a typical Stark map. As the electric field is increased, 
different energy levels are shifted by different amounts. Adiabatic (dotted arrow) 
and diabatic (dashed arrow) paths are shown. The inset shows the structure of an 
avoided crossing, with energy gap E a . 

this energy gap, it is likely that the avoided crossings are traversed diabatically 
(dashed arrow in fig. [1]), and the atom is ionized at a point at or above the 
classical ionization limit. Intermediate slew rates, however, result in a random 
combination of adiabatic and diabatic traversals, leading to the ionization 
field becoming much less well denned, which significantly reduces the state 
discrimination of the detector. The states of interest in our experiment, around 
n = 60, are in this intermediate regime and therefore effective detection has 
been a persistent problem. 

The Rydberg states have an average lifetime of the order of 0.1-lms, so it 
is also necessary to complete the field ramp before the atom decays from 
the Rydberg state. The detector must also fit inside the the physical space 
available inside the experimental vacuum chamber. This places constraints on 
both the maximum and minimum timescales for the different measurements. 

Traditional designs of field ramp for single state selective detectors have tended 
to involve an approximately linear ramp, with more complicated combinations 
of linear ramps of differing slew rate to distinguish between different states. 



4 Technical details of the experiment 

The experimental arrangement controlled by the GA is shown in fig. [2J A 
beam of rubidium atoms are excited to high-lying Rydberg states via a three 
step laser excitation process. These atoms pass through an ionization region 
where the ionization is driven by the electric field between a pair of parallel 
plates. The liberated electrons are collected in a channel electron multiplier 
(CEM) which records a TTL pulse for each detection event. The atomic beam 
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Fig. 2. The experimental setup that was optimized by the GA. The laser excitation 
is a sequence of three infrared lasers at 780nm, 776nm and 1257nm. A time varying 
voltage is applied to the field plates and ionized electrons are collected at the channel 
electron muliplier (CEM). 

and detectors are contained within a vacuum chamber at a pressure of around 
10 _8 mbar. The laser light is coupled into the chamber via optical fibers. 

This arrangement is a subset of the quantum optical micromaser experiment 
[8] that probes fundamental quantum mechanical effects, such as the con- 
trolled creation of arbitrary numbers of photons [H] and entanglement |l(Jj . 
The frequencies of the first two laser steps are locked to the relevant atomic 
transitions inside a vapor cell on the optical table. It is not possible, however, 
to lock the third step in this way, and for this experiment this laser was locked 
with an optical frequency synthesizer (Menlo FC1500). 

4-1 Temporally varying vs. spatially varying profiles 

For this particular experiment, we chose to use a temporally varying field 
(sometimes referred to as pulsed field ionization [7]) since it is technically eas- 
ier to achieve, requiring only a single electrical feedthrough into the vacuum 
system. To ensure that we are only detecting atoms that are in the correct 
portion of the field ramp, it is necessary to gate the electron multiplier's out- 
put. Clearly this duty cycle puts a limit upon the maximum count rates that 
will be measured. 

4-2 The Problem to be solved 

The first experiment that has been performed, as a proof of principle, involved 
simply a beam of atoms prepared in a particular Rydberg state (in this case 
a state with n = 63) impinging upon an externally controlled field ionization 
region with a single electron multiplier. 
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Fig. 3. Pulsed excitation scheme. The laser excitation is pulsed on for a short period, 
exciting a small group of atoms moving at a wide range of different velocities. Only 
those atoms that are ionized during the gate pulse are recorded, hence selecting 
only a very narrow velocity class. 

Broadly speaking, there are two different scenarios that a detector will be 
required to cope with. First of all, many cavity QED experiments require a 
fixed and very precise atomic velocity, requiring a detector to be optimized 
very tightly to that velocity, perhaps at the expense of other velocities [TT] . 
The majority of experiments will, however, require the ability to detect atoms 
moving at a wide range of velocities. For example, a typical experiment will 
consist of a sequence of atoms travelling at identical velocities, followed by 
one or more travelling at a different velocity as a probe of the quantum state 
inside the cavity [9]. 

To optimize for the first type of experiment, we have two schemes for velocity 
control at our disposal. First of all, we can align the laser excitation at an angle 
to the atomic beam to ensure we only excite atoms with the correct Doppler 
shift (and hence correct velocity). The second scheme is a time of flight style 
system where the laser excitation is pulsed and the CEM gate occurs at a 
fixed time afterwards. This type of scheme is illustrated in fig. [3l 

For the second type of experiment, we require the evolved solution to be ro- 
bust to varying atomic velocities. This is achieved by simply preparing the 
full Maxwell-Boltzmann distribution of atoms and optimizing the detector si- 
multaneously to all velocity classes present in the atomic beam. To do this 
we apply a continuous perpendicular excitation to the atomic beam and gate 
the electron multiplier synchronously with the evolved waveform, as shown in 
fig. HI This allows us to gate the detector many more times a second and eval- 
uate the fitness of a given individual with correspondingly higher confidence. 
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Fig. 4. Continuous excitation scheme. The laser excitation is applied continuously 
so that atoms of all velocities are present in the detector during the gate pulses. 
This scheme allows much faster repetition rates and hence offers better counting 
statistics. 

A possible combination of the two schemes would be to create a catalogue of 
profiles optimized for different velocities and actively switch between them for 
the velocity of the next atom. 

A spatially varying profile has a number of advantages over a simple tem- 
porally varying profile. The main advantage is that there is no need to gate 
the electron multiplier's output signal, since every atom experiences the same 
field, regardless of the time at which it arrives. This immediately gives us a 
gain over the limited duty cycle necessary in a temporally varying system. 

There is also a subtle difference in that the spatially varying system will pro- 
duce a slew rate that is a function of the atomic velocity, unlike the temporally 
varying system, where a given instantaneous electric field value is experienced 
at all positions inside the detector. 





5 The Genetic Algorithm 



The genetic algorithm used in this case is based upon the Microbial GA |12j . 
In summary, this scheme of GA utilizes a pairwise tournament selection, the 
parametric uniform crossover is implemented by overwriting elements of the 
tournament loser with the corresponding elements from the tournament win- 
ner, with probability P c = 0.5. A standard creep style mutation is employed, 
with the mutation size X being drawn from a normal distribution centered on 
(X) = 0V with standard deviation ax = 5V. 
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The initial population consisted of 50 individuals of 20 evolvable vertices that 
were initialized with random values uniformly distributed over the range of 
-45V to 25V. 

5.1 Fitness Function 

A key element of any GA is the choice of fitness function that drives the selec- 
tion procedure. In this simple example of the experiment, the fitness function 
is chosen to be the absolute count rate, since a higher count rate indicates 
a higher detection efficiency. Future designs will simultaneously optimize for 
two different atomic states, and hence will require a correspondingly more 
complex fitness function. A key advantage of a tournament style of GA is that 
the pairwise selection helps to negate any effects of medium to long term vari- 
ations in the experiment. For example, slight laser drift may cause the flux of 
atoms in the 63P state to fluctuate, which would make a generational scheme 
unreliable, as the best individuals won't necessarily achieve the best fitness. 
Since these types of variations typically occur on a timescale longer than the 
time to evaluate just two individuals (a few seconds), we can be confident that 
the tournament winner actually is the better individual. 

5.2 Genotype 

The genotype that was evolved was a set of vertices that defined a waveform 
which was fed to an arbitrary waveform generator {Tektronix AFG3022). The 
voltage output of this device was then amplified and fed to the parallel elec- 
trodes inside the vacuum chamber (fig. [2]). The time varying voltage applied 
to these plates defines the electric field profile in the detection region. 



6 Results and Analysis 

The GA was allowed to run for approximately 2000 tournaments, with each 
evaluated individual being recorded along with its fitness. The results are 
plotted in fig. [51 showing that the expected fast field ramp at the time of 
the detector gate (75-85us) appears very early in the evolution, as might 
be expected. More unexpectedly, further ramps appeared in the field profile 
before the detector gate. The presence of a single negative dip before a large 
positive pulse has been noted before in a similar experimental system |13j . 
the explanation for which is that the atoms are manipulated into a state that 
preferentially ionizes exactly at the classical limit. The presence of multiple 
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Fig. 5. The winner of each tournament during the evolution. The initial random 
population quickly develops a series of sharp oscillations followed by a plateau where 
the detector gate is active, between 75-85us. 

dips before the main pulse in our experiment suggests that the system is 
seeking a state which is more easily ionized. The more complicated structure 
most probably occurs due to the different value of n we use, placing us in 
the intermediate slew rate regime as opposed to the much higher n used in 
reference [13J, which is firmly in the purely diabatic regime. 



7 Future Directions 

The first extension is to incorporate a second detector and optimize for state 
selectivity. This would involve an extended fitness function that provides se- 
lection pressure to increase the efficiency of each detector individually and also 
for their ability to distinguish the correct state. The system must be cooled 
to cryogenic temperatures (as it is in the micromaser experiment) in order 
to suppress thermally induced intermixing of the Rydberg states that would 
introduce noise into the fitness evaluation. 

In typical experiments, the experiment is run continuously, so a pulsed detec- 
tion scheme is unsuitable. In this spatially varying field is required. To 
accomplish this, we are in the process of building a new detector where the sin- 
gle pair of electrodes are replaced by an array of strip electrodes. The voltages 
sent to each of these electrodes will then be the genotype to be evolved. 
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